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A slave boson representation for the degenerate Hubbard model is introduced. The location of the 
metal to insulator transition that occurs at commensurate densities is shown to depend weakly on 
the band degeneracy M. The relative weights of the Hubbard sub-bands depend strongly on M, as 
well as the magnetic properties. It is also shown that a sizable Hund's rule coupling is required in 
order to have a ferromagnetic instability appearing. The metal to insulator transition driven by an 
increase in temperature is a strong function of it. 
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There has been dramatic progress in our understand- 
ing of the Mott transition in the last few years. Care- 
ful experimental studies of systems iruthe vicinity of the 
Mott transition have been carried outEl and two new the- 
oretical Jjools, slave bosons mean field theories (see for 
instanceB and references therein), and the limit of in- 
finite dimensions have been adapted to its study. For a 
review sea3. Most of the modern work has focused on the 
single band Hubbard model. Now that both the metallic 
and the Mott insulating states of the (doped) titanates 
and vanadiates have been studied experimentallyuu (cor- 
responding to 3dl and 3<i2 configurations in the Mott 
insulating state), there is a need for a theoretical frame- 
work allowing for understanding the Mott transition for 
arbitrary degeneracy and density. This paper is aimed 
to provide such a technique and to apply it to a variety 
of quantities that cannot be obtained easily using alter- 
native approaches. Most of the results are obtained in a 
closed analytical form, allowing for a qualitative under- 
standing of the physical situation. 

In this work we investigate the effect of strong Coulomb 
interaction in systems with orbital degeneracy. Such a 
situation is realized in virtually all transition metals and 
transition metal oxides. These systems contain d elec- 
trons in cubic or trigonal environments, the crystal field 
can only lift partially the degeneracy of the d-bands, 
down to 2 as is the case of V2O3EI or 3 as in LaTiOz. 
Our goal is to understand how degeneracy affects the be- 
havior of the different physical quantities near the Mott 
transition. To carry out the investigation we extend the 
slave boson technique which has been very successful in 
the study of the Mott transition, to the orbitally degen- 
erate case. Compared to the variational wave function 



approachta our formalism is more flexible since, as we 
demonstrate in this paper, it allows us to calculate a 
variety of quantities which are not easily accessible to 
the variational approach, as a function of the correlation 
strength and doping. It can also be improved systemati- 
cally by performing a loop expansion around the saddle 
point. Our main results are the following :a) low energy 
single particle quantities such as the critical value of the 
interaction strength at which the transition occurs, the 
quasiparticle residue and the single particle Mott Hub- 
bard gap depend very weekly on degeneracy justifying 
the agreement between theory and experiment when it 
was applied to orbitally degenerate systems, b) the rel- 
ative weights of the Hubbard bands depend strongly on 
degeneracy in agreement with other methods. El c) the de- 
generacy temperature decreases with increasing band de- 
generacy d) the magnetic properties, the magnetic sus- 
ceptibility and its associated Landau parameter in the 
paramagnetic phase and the magnetic phase diagram is 
strongly modified from the one band case. 

The Hamiltonian describing the low-energy properties 
of these systems is commonly written as: 

H = ^tp^^o + ^3 n hP,1 n i,P,l 

i,j,a,p i,p 

i,p'=£p i,a,p'<p 

where a is a spin index for the up and down states while 
p is labeling the M bands, and U n = U + nJ. Taking J 
finite accounts for the Hund's rule coupling. 

As for any model with on-site interaction, a slave bo- 
son representation can be introduced, mapping all the de- 
grees of freedom onto bosons. We can re- write any atomic 
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state with the help of a set of pseudo-fermions {/«} and 
slave bosons {ipau...a m } (0 < m < 2M). ^},„a m 

is the 

slave boson associated with the atomic state consisting of 
m electrons in states \a±, a m > where a is a compos- 
ite spin and band index. By construction it is symmetric 
under any permutation of 2 indices, and if any 2 indices 
are equal. We can now write the creation operator of a 
physical electron in terms of the slave particles as: 



c+ = z+f+ 



2M 



i J a,w]...,am-i^arL,...]a m --L ai ^ 



a 



m—l ai<.<a m -i 



Ro 



2A/-1 

E E 

m— Qi,<. <,atrt 
2M 



/-„ [i E ^ 



m=l ol\ < . <a n 



Yoti,.,a m 'f J ai,.,a m l u 'r 
^(m) l-i 



a 



Namely L a normalizes to 1 the probability that no elec- 
tron in state \a > is present on a site before one such 
electron hops to that particular site, and R a makes sure 
that it happened. Clearly the eigenvalues of the opera- 
tors L a and R a are 1 in the physical subspace. Now, the 
redundant degrees of freedom are projected out with the 
constraints: 



We now proceed to the mean-field theory, and we inves- 
tigate the paramagnetic, paraorbital saddle-point. The 
latter is obtained after integrating out the fermions, and 
setting all bosonic fields to their classical value. The Mott 
transition that occurs at commensurate density n is best 
discussed by projecting out occupancies that are larger 
than n + 1 and smaller than n — 1 (if any). The con- 
straints allows for eliminating the variables tp( n ~ 1 > and 
ipW to obtain the grand potential at n: 



(2) 



Q{D) = (1 - 2D 2 )D 2 {^Kli + V^feo 



z+ describes the change in the boson occupation numbers 
when an electron in state a is created as: 



U(D 2 



) - MP 



(7) 



2M J deep(e)f F (z 2 e + A - /i), D 2 



and 



(3) 



The operators z+ in Eq. ([?]) describes the change in 
the slave boson occupation as a many channel process. 
In order to recover the correct non-interacting limit at 
mean-field level, one has to observe that the classical 
probability for these processes to happen is not simply 
given by taking the Bose fields in (H) to be given by their 
classical vakies, but by introducing normalization factors 
L a and i?J§0 as z+ = Y^'<l> +(m) L a R a ip {m - 1 \ where: 



with 6o 

(2M^(„+l) 2) M = (2M _ n + 1)/(2M _ ^ 

c n = (n + l)/n. Minimizing Eq. (0) with respect to D 
yields a critical interaction strength at which D vanishes. 

r( "' M) - -eoiy/b— 1 ^ 2 



It reads Uc = — eo( \J b n ,M + \J~Cn) which reproduces 
the results of the Gutzwiller approximationEM This lo- 
cates the Mott transition. Restricting ourselves to a flat 
density of states we can relate the critical interaction 
strength to the band width W . We obtain: 



jj(n,M) 



nW 
4M 



(2M -n)(^/b nM + Vc^) 



(8) 



Its band degeneracy dependence is fairly weakfcl. The ef- 
fective mass of the quasi-particles diverges at the Mott 
transition. We obtain: 



nr 



r (n,M)2 



(9) 



8 m 

^ Due to the particular form of the coefficients b and c the 
dependence on the band degeneracy is weak. The crit- 
ical interaction strength increases with M so the quasi- 
particle residue Z increases slightly with M. For small 
values of U, (which we treated without projecting out 
higher occupancies), Z decreases with increasing M. So 
there is a crossover value of the interaction strength be- 
yond whichpthe system becomes more metallic with in- 
creasing MxQ. As a function of the hole doping S, the 
quasi-particle residue vanishes for 6 going to above 



f+f— X^ 2M \" ilM m ) i/; (m) = 

J a J a L-im=\ / / ra,a 1 ,.,a ra _ 1 ra,ai,.,a m _i u 

ai<.<tt m -l 

EZo E ^L^^-l-O . (5) 
£>l<.<a m 

We obtain the Lagrangian at J = as: 

ia i,m ia.i<.<ot m 

/ \ m 

+ ^iJ Z t,afi^a Z j,otfj,m,a ■ (6) 



U^ M) as: 



Z = ~{K,M - C n ) + -Z-((bn,M + C n ) 



X y/l + ±ip n ,M + 4 v /^n,MC„^„,j\/) 



where we introduced: 



U, 



(n,M)2 



(u-uh M) )(u-m ,hxvi '{ 



) 2 ) 



(10) 



(11) 



The expression of the quasi-particle residue consists of 2 
contributions which arc cither symmetric or antisymmet- 
ric with respect to particle or hole doping. The antisym- 
metric contribution vanishes for n — M as a consequence 
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of the particle-hole symmetry. The asymmetry of z 2 on 
particle or hole doping is seen to increase under an in- 
crease of | n — M | . It vanishes more slowly for hole doping 
(for n < M) than for particle doping, for increasing de- 
generacy at fixed n, for increasing degeneracy at n = M 
and under an increase of U. As an example we calculate 
the effective mass for the 2 band model and show it on 
Fig. [j], which has been calculated without projecting out 
higher occupancies. 

Effective Mass 

M=2 



U/U c {2,1): 
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FIG. 2. Chemical potential for n = 1 for the 1 band 
(dashed line) and 2 band (full line) models. 



(n,M) 



the Mott gap 



(n,M) 



as A 



In the limit of U » U, 
is given by U, while it closes at U, 
jj{n,M) I jj/jj( n ' M ) _ ^ the square root behavior being 
typical of slave boson mean-field theories. It can be read 
from Fig. |2| where itJs compared to the 1 band case as 
obtained by LavagnaHl. Clearly going from 1 band to 2 
bands does not imply a big difference in the Mott gap. 
Indeed we obtain the A/t r i™' M ' ) is independent of M at 
n = M, while for fixed n it depends very weakly on M. 

Our result can be compared to experimental data. For 
the series La x Y\^ x TiO^^ Okimoto et aO measured how 
the gap depends on the band width. Assuming (for large 
ratio U/W) A ~ U — W we obtain out of their data 
U = 3.2eV. Inserting this and the experimental value 
of (U/W) c ~ 1.3 in Eq.(|l|) we can compute A/W as a 
function of W/U and compare it with experiment on fig. 
|j| The experimental trend is clearly reproduced and the 
quantitative agreement is very satisfactory. 
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FIG. 1. Inverse effective mass in the 2 band model as a 
function of density for several values of U. 

Interestingly we also obtain a Mott gap. Indeed the 
number of quasi-particles is a continuous function of their 
chemical potential p — Ao/2. However the saddle-point 
equations show that the Lagrange multiplier A jumps 
when going through the Mott gap which implies that A 
is going to jump as well, and so does p. As a result we 
obtain the Mott gap A = lim^Q- n(5) — lim^^o+ /i(<5) 
as: 



A = J(U - U^ M) ){U - ut M \^^ 



:) 2 ) 




(12) 



FIG. 3. Dependence of the Mott gap on the band width for 
n = 1 and M = 3. Circles: experimental data of Okimoto et 
aP. 
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We now turn to the Hund's rule coupling dependence 
and treat as an example the two band model around the 
n = 1 Mott insulating lobe. At p = 1 the grand potential 
at the saddle-point reads: 



n = -e (1 - 2r 2 ) (r + (do + d x + A )/\/2) 2 
+ (U + 3 J)A 2 Q + {U + J)d 2 x + Udl - pp . 



(13) 



with do = 
A = 



(v>So 



(2) 

tit 



+■ v<S)/V2, 



d x = 
: = dl 



A 2 



and 



A = A Q /2, and we have used the constraints to re- 
move the variables an d ip^ 1 '. 

Such an expression differs from an ordinary Ginzburg- 
Landau free energy in that respect that it cannot be writ- 
ten as a 4th order polynomial in the variables do, d x and 
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Aq. As a result, if there were to be a critical point for 1 
field, it would be critical for the other ones as well. We 
obtain the location of the Mott transition as: 

^ = 1 -| ^+°( J / C/ ) 2 ) • ( l4 ) 

Another regime of interest is the large J regime. There 
we obtain the location of the Mott transition as: 

= - jje (3 + 2V2)(1 - I'j) + 0{Cjf) (15) 

and thus decreasing J from oo leads to an increase of the 
critical interaction. Another intriguing feature of transi- 
tion metal oxides such as V2O3, is the metal to insulator 
transition that occurs in the vicinity of the tri-critical 
point under an increase of temperature. It has recently 
been interpretatedtj as the transition from a Fermi liquid 
with finite quasi-particle residue Z to an insulator with 
Z = 0. In other words there is a finite coherence tem- 
perature T co h at which the coherence of the Fermi liquid 
(and Z) vanishes. This result was obtained in the dy- 
namical mean-field approximation to the 1 band model, 
which becomes exact in the limit of large dimensions, and 
recovered in the Gutzwiller ApproximationE-1 At finite 
T there is a first order metal to insulator transition at a 
critical U { C M) {T): 

U { C M) (T) = U { C M) (0) - \lw { c M) {0)T\n2M (16) 

Thus an increase in temperature may produce a metal to 
insulator transition, which is consistent with the exper- 
imental situation in V2O3. In the dynamical mean-field 
approximation at finite temperatures there is an interac- 
tion strength U C 2(T) at which the metallic solution ceases 
to exist. This quantity can also be evaluated in our slave 
boson scheme and is given by: 

ug°(T) = C/i M >(0)(l - a M (T/W)i) (17) 

with ot\ ~ 2.53 and ai ~ 3.32. 

We now turn to the calculation of the magnetic sus- 
ceptibility. Here we generalize the calculation of Li 
et a£2l to the 2 band model. The linear response to 
an external magnetic field is obtained as a 1 loop cal- 
culation of the correlation function of the slave boson 
fields in the spin-antisymmetric band-symmetric chan- 
nel. 3 fields couple in this channel!!]: %_ = | J2 a a,l Pa\ 

X+ = T^t!? _ ^u) ^ K = |Ea ^a. and the 

magnetization is expressed in terms of slave bosons as 
M = 4rfox++2?iU' ) X-- The resulting susceptibility arises 
as an RPA formliH 

XS = l + F$xo/N(E F ) ■ (18) 

We now determine the instability line of the paramag- 
netic phase with respect to ferromagnetism. For J = 
we find no ferromagnetic instability even near the Mott 



transition, while for finite J we find that the Mott metal 
insulator transition may be preempted by the appearance 
of a ferromagnetic phase. In other words, a sufficiently 
strong Hund's rule coupling turns a Mott insulator into 
a ferromagnet. Originally the Hubbard model was in- 
troduced in order to describe ferromagnetism in narrow 
band systems, but it has been recently established that 
the ground state is not ferromagnetic for apy. reasonable 
values of the parameters on generic latticeslIjO. We find 
that the ground state is much more likely to be ferro- 
magnetic in the degenerate model for finite J, as shown 
in Fig. § 
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FIG. 4. Instability line of the paramagnetic phase for 
U/J = 10 (dashed line) and U/J = 5 (solid line). The di- 
amond (square) indicates the position of the Mott transition 
for U/J = 10 (U/J = 5). 

Our method can be applied to the calculation of dy- 
namical quantities too. In the strong coupling regime 
the one-particle excitation spectrum is split off into sev- 
eral pieces, each of them carrying some fraction of the 
spectral weight (which are adding up to 1 so as to ful- 
fill the sum rule). The various pieces are following from 
the discrete atomic levels, which are well separated by 
multiples of U, broadened by exchange processes. Let 
us now determine the fraction of the spectral weight car- 
ried by each sub-band. In our language the low energy 
excitations are involving the field and the high en- 
ergy excitations centered around U the field ip^ 2 \ Higher 
energy excitations involving higher local occupancies are 
left out. We obtain the spectral weights in both bands 
from the decomposition of the physical electron opera- 
tor (Eqs. (§,§))■ Accordingly the spectral weight of the 
Green's function T < c a (r)c+ (0) > in the lower Hubbard 
band (Wlhb) an d the upper Hubbard band (Wjjhb) are 
given by: 

and are shown in fig. 0. Here the weights do not quite 
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add up to 1 in the 2 band case because we projected out 
occupations larger than 2. In other words, on top of the 
2 sub-bands which are considered here, there appear a 
second upper Hubbard band (centered around 3U — 3/i, 
corresponding to triple occupancy) which is becoming 
relevant in the particle doped regime. To a very good 
accuracy its contribution to the spectral weight is given 
by 1 — Wlhb — Wuhb- Clearly the degeneracy plays an 
important role as the weight of the upper band at n = 1 
in the strong coupling regime is given by (2M — 1)/2M. 
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FIG. 5. Spectral weight of the upper (dashed lines) and 
lower (solid lines) Hubbard bands for the 1 band and 2 band 
models, at U = 2U, 
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In summary we introduced a slave boson representa- 
tion of the degenerate Hubbard model. We obtained that 
the band degeneracy has a weak influence on the location 
of the Mott transition, while the degeneracy temperature 
and the dynamical and magnetic properties strongly de- 
pend on it. We also showed that no ferromagnetic in- 
stability occurs unless the Hund's rule coupling becomes 
sizable, yielding a generic scenario for ferromagnetism in 
transition metals and transition metal oxides. In that 
case a ferromagnetic instability may even shadow the 
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